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The electrostatics of 2D system with complicated inner boundary is studied. The object which we 
call "monster" is built by an iterative process of multiple conformal mapping of the circle exterior. 
The procedure leads to the figures built from circle with branching curved cuts (i), to the multiple 
tangential near-round circles (ii) and the Mandelbrot map (iii). The polarizability of a monster in 
the homogeneous external field is found. The light scattering cross-section was expressed through 
the polarizability of a monster. 

PACS numbers: 78.20.Bh, 78.66.Sq, 05.45.Df 

Introduction. The fractal aggregates appear in many experimental and theoretical situations. Such aggregates 
are the topic of the diffusion limited aggregation, the theory of percolation and the electric breakdown. 

The purpose of the present paper is to study the polarizability of solitary particles with complicated boundaries. 
We shall refer the complicated particles, built from small subunits, as "monsters", distinguishing them from fractals, 
which result from the infinite growth. The monsters can be obtained both directly in the growth process, and as a 
result of Van der Waals aggregation of preproduced spherical nanoclusters into large snowballs. It is natural to expect 
that in many real situations medium is composed from the monstrous or fractal objects, rarely embedded into matrix. 

The optical properties of fractals, including polarizability were studied in Ref. pH|. The attention of these studies 
was mostly concentrated on the localization of high-frequency polar cigcnmodes of fractals and corresponding large 
^ ' fluctuations of local electric fields. 

The works,Di3'Q most close to the purpose of the present study, consider the fractal as a molecule built from counled 
dipolar monomers. The system of linear equations for dipoles was solved either in self-consistent approximatiorJj or 
by computer simulations.l3tl 

i^i - The self-consistent approximation is valid only if the polarizability of the fractal medium is not large. The large 
polarizability of the medium leads to the strong screening of the external field. This is the case for metal clusters, 

■ studed in the present work. 

Unlike the referred papers, we are interested in the situation where the local fields are strongly deviate from the 
mean value. We shall consider the aggregate interior as a macroscopic medium, neglecting the atomic structure effects. 
In this approach one should solve the macroscopic Maxwell equation inside and outside of fractal aggregate with some 
boundary conditions on its surface. We shall deal with aggregates, that are small in comparison with the wavelength. 

■ This permits to solve the Laplace equation, instead of the Maxwell one. 

Q\ \ The Laplace problem is significally simplified, if one does not need the solution both inside and outside of the 
0^ • object. It is so, if the dielectric permittivity of the cluster is essentially larger or essentially smaller than that of the 
external medium. If the cluster medium is metal, then its permittivity is determined by the formula 1 — uji/uj 2 , where 



pi 

u) is the light frequency, and ui p is the plasma frequency. The first case corresponds to ui <C u> p , the second one to the 
frequency near the bulk plasma frequency. 

We shall study the 2D case, which allows to use the theory of analytical functions for solution of the Laplace 
equation. This simplification may be considered as a model for the 3D case. At the same time, it is applicable to the 
q | clusters with cylindrical shape F(x, y) = const. 

We shall find the electric polarizability of a large aggregate of complicated shape that can be constructed by means 
of multiple conformal maps. In this case the exact solution of the problem can be found. One kind of such construction 
of 2D fractal objects was suggested in Ref. [To|-|T^. The fractal was treated as a result of iterative mapping of a simple 
domain, for example, exterior of a circle. The elementary step maps the exterior of the circle onto the exterior of a 
circle with a bump or with a strike. The multiple mapping complicates the form of the domain, producing a bump 
(strike) for a step of an iteration. This construction permits to simulate growth, guided by so called harmonic measure, 
when the addition of an element of fixed area to the domain appears with the probability, proportional to the diffusion 
flow, determined by the solution of the harmonic equation of the diffusion (DLA). ._. 

Another approach to constructing fractal is a simple map, for example, square mapuJ z n+1 = c — z^. The multiple 
superposition of this map separates the complex plane into the repulsion basin, composed from all points zq for which 
z n — > oo. The boundary of that basin is a fractal, called Julia set. 

The limited repetition of maping produces a monster, the complicated but not fractal object. We are interested in 
the properties of the monster as well as the properties of the limit of the infinite number of repetitions, the fractal. 
The complexity of the fractal grows linear with the map order. It differs this construction from the square map of 
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Mandelbrot, where the complexity grows exponentially with the map order. 

In the case of rare inclusions the effective electrical properties of medium are determined by individual monster 
properties. For example, the effective dielectric permittivity of the medium, containing monsters with concentration 
n and polarizability Xij m the limit of low concentration is 



(1) 



Similar formula cry = a(5ij + Aimxij) determines the effective conductivity of the medium with conductivity a, 
containing insulation inclusions. 

The other important property of physical structures is the light scattering cross-section. It is determined by the 
polarizability, if the wavelength exceeds the diameter of the object. In particular, we shall deal with the polarizability 
and with the scattering of light on a small metal particle. 

Method. These properties can be found by solving the Laplace equation <f> A<fi(r) = for the potential outside 
of a solitary particle in the external uniform electric field E = E x + iE y . For lu <^ oj p the boundary conditions at the 
particle surface S and at the infinity are 



d S (b = 0, 



-Re{E*w). 



(2) 



Here w — x + iy. 

The effective conductivity of media with rare dielectric inclusions is determined by the solution of the current flow 
problem around the inclusion. 

If oj — ui p <C lOp, the potential ip satisfies the condition d n ip\s = at the boundary of the inclusion. According 
to the Couchi-Ricman conditions <p and ip arc the real and imaginary parts of the complex potential <I>. This case 
corresponds to the problem of the current flow in a conducting medium around an insulating inclusion with the same 
shape. 

The idea of the present work is to construct a complicated object from simple one by means of repetition of some 
conformal map f(z). The map should transform the exterior of the basic object to the exterior of complicated one 
and does not move the infinity point. It doesn't alter the Laplace equation and the surface boundary conditions (with 
the accuracy to the coefficients). Therefore, the solution of the problem for the source object gives it for the image 
object too. 

Let us consider a multiple map f( n ) composed of the n basic maps f\(z), where A is a parameter (or a set of 
parameters) of the map, and some additional map G: 



F(">(™) = G°/ Al °/ A2 °...°/ A „H. 



(3) 



The parameters Xi may be chosen fixed, random, or depended regulary on i. 

If we are interested in the dipole component of electric field only, we should take into account only /i and /_i 

coefficients of series expansion of the map f\i(w) = fi w + + /i^w -1 + .... 

The complex potential outside the unite circle in the external uniform electric field in 2D case is 



$H = -(E*w- E/w), 



(4) 



on the complex w-plane. Let G = 1. Substitution of the series expansion of the inverse function w(z) : F^ n \w(z)) = z 
into (H) gives us the polarizability tensor x- 



X= 2 



A ( Re{B) - a Im(B) 

Im(B) —Re(B) - a 



(5) 



where a = — 1 in the case of metal particle and 1 in the case of dieletric one. A and B are defined as follows: 
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B =A 



fe-i 



k=2 



i=l 



(6) 
(7) 



The principal values of polarizability tensor are 



Xl,2 



~{°±\B\) 
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Examples. In the recent worksliH~£3 a multiple conformal map (i) was used for constructing of models of fractal 
objects, in particular, DLA and dielectric breakdown clusters. The map (@) transforms the unit circle to the unite 
circle with a small bump (or strike) on it (Fig. fy). 



f x ( w )= w ^ll±±(l + w ) 



1 + 10 + 



1 2 1 - A 
+w | 1 + — — — - 

w 1 + A 



w 



1/2 



(8) 
(9) 

The shape of the bump is determined by two parameters a and A. In this work we assume them to be constant. We 
use the rotated map (^|) as the basic for F and consider simple dependence of the angle of rotation 6k — k9, where 
9 = const is the angle incommensurable to tt. This choice enables us to find the exact equation for the polarizability 
but still produces the complicated objects. The map, recursively applied to the source unit circle, produces monsters 
as on the Fig. [[} 

The coefficients A, B of the polarizability tensor in the limit n — > oo are: 

A = (1 + A) 2a,i , 



B = aX 



2 + (2a- 1)A 



(1 + A) 2 P 2 - 2/3 cos 26» + 1 



(10) 



where /?=(! + A) 2 




FIG. 1. The case (i). The images of the unit circle according to the Eq. g after n = 50 iterations for parameters 9 — 2.7, 
a — 2/3, A = 0.1 (case of strikes, figure a) and 9 = 2.15, a = 1, A = 0.1 (case of bumps, figure b). Images of basic maps applied 
to the unit circle are shown in the corners (the size of basic map is arbitrary) . 



Another considered basic map (ii) is the map of the exterior of the unit circle to the exterior of two tangenting unit 
circles: 



fW 



log 



w — 1 



+ 1, f e (w) = e i6 f(e- ie w). 



(11) 



The example of multiple map is shown on Fig. || in the case of 9^ = k9. 
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FIG. 2. The case (ii). The image of the circle after 10-th map ( |ll"| ) for 9 n — n8, 9 — —2.1. The basic map is shown also. 
In the case of n — > oo corresponding quantities A and B for polarizability (|^) are 



2n 



£ = 



7r 2 (-e 2 ^) 



3(4 + 7r 2 e 2ie ) ■ 

The third basic map (iii) is the known Mandelbrot map. Resulting transform is 

F("'H = (1 - (1 - (... - w 2 ) 2 ...) 2 ) 2 '" . 



(12) 



(13) 



The additional transform G(a) = a 2 makes the final map single-valued on the exterior. On Fig. the domain is 
shown, which boundary, so-called Julia set, is the image of a unit circle. The multiple map (iii) produces fractal, which 
appears, for,example, in the hierarchical description of the percolation transition in the complex plane of percolation 
probability.^ 



FIG. 3. The case (iii). The Mandelbrot fractal, obtained by the map (jl3|) of the unite circle. 
The polarizability of the metallic Mandelbrot monster is 

/ r- \ 2 "" +1 

(1 + V5)\ 



Xxx,yy — i 



! Xxy,yx 0. 



(14) 



Light scattering. The 2D problem may be considered as a model for the more realistic 3D system. From the other 
hand, it describes the properties of the cylindrical objects, which size in z-direction exceeds the maximal diameter in 
(x,y) plane. 

In relation to the problem of light scattering we shall consider the geometry, depicted on the Fig. ^. In this case, if 
the electric vector E = |E|e lyes in the face plane of the monster, the cross-section is 



^,4, |2 4 sin (¥0 



(15) 



where is the vector of polarization of the electric field. 



} 



FIG. 4. The geometry of the light scattering problem on a cylindrical object with complicated section, k is the wave vector. 
The plane of polarization coincide with the cross-section plane. 
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Discussion. The polarizability of the round metal particle in 2D is roughly proportional to the second power of 
its radius. For a strongly extended asymmetric particle the components of polarizability tensor are determined by 
the sizes of the particle in the direction of the field. 

In general, the monsters in the cases (i) and (ii) have no symmetry so each component of tensor x is finite. In the 
case (iii) the fractal has x and y axes of symmetry resulting in Xxy — 0. The monsters (i) and (ii) grow exponentially 
with the number of iteration. This leads to the same growth of the polarizability. 

In the case (iii) the sizes of monster are limited, as is the polarizability tensor. The limiting ratio of its components 
roughly corresponds to the ratio of sizes. The fractal structure of this system manifests itself in the exponential law 
of approaching to the limit. 

The purpose of the present work was to find the polarizability of metal clusters. However, these results are more 
widely applicable. They describe the effective conductivity of a medium with insulating or perfectly conducting 
inclusions. In accord with the Einstein relation, the formulae similar to (|l|) determine the effective diffusion coefficient 
of some particles in the medium with inclusions, non-penetrable for diffusing particles. 

We have used some ways of constructing monsters and fractals by means of multiple conformal maps. The basic 
map may be more simple as the Mandelbrote square map (iii) or more complicated as the bump map (i) or the 
tangental circle map (ii). The advantage of complicated maps (i) and (ii) is that they permit to adjoin unites of fixed 
form to a cluster, the disadvantage is that the complexity of monster grows linearily with the iteration number, while 
the complexity of square map (iii) (and more general one, = c — 2„) grows exponentially. 

We have considered the sequent maps with fixed parameters only. In the cases (i) and (ii) the external parts of 
monsters are essentially larger than internal ones. The fine details of the monster are inside of it and the exterior of 
the monster is not fractal. The polarizability of the monster depends generally on its external shape, scaling like the 
square of its diameter. This differs situation from one, considered in Ref. [10-12, where the fractal shape of a growing 
cluster is achieved by the selection of the map parameters Xi in such a way that adjoining bumps have equal sizes. 
Our results may be easily generalized to this case also and to the other cases of random, or some regular dependence 
of parameters on the map order. 

The specific shape of monsters (i) is similar to the branching traces of electric breakdown and hence is applicable 
for description of electric properties of medium after electric breakdown. The monsters (ii) have the shape of coupled 
droplets in emulsion. By selection of parameters this case may be generalized to desribe clusters of near round droplets 
with equal sizes. 
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